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The ID phase diagram of a model for correlated hopping of electrons in a lattice of Berry phase 
molecules is presented. Electrons hop in presence of an extra orbital degree of freedom at each 
site. This is mimicked as a spin-1 variable whose allowed states depend on the electron occupancy 
so as to take into account the orbital degeneracies of different molecular occupancies. In the ID 
case we find that at low electron densities n << 1 there is a region with dominant superconducting 
correlations surviving an additional repulsive on-site interaction U as strong as the bandwidth, 
W = At. The critical value Uc of U below which superconductivity is found to be dominant decreases 
with increasing density n. For n = 1/2 we find Uc/t « 1, whereas at n = 1 (half-filling) our (less 
accurate) results are compatible with Uc/t ~ 0. For U > Uc{n) and away from half-filling (n 7^ 1) 
the system is metallic with dominant 2fc_F charge density wave (CDW) correlations. At half-filling 
a charge gap opens for U > Uc and the system becomes an insulator. A spin-gap characterizes the 
phase-diagram for all densities and for all values of U, even in the metallic regime U > Uc. 
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I. INTRODUCTION 



In the normal theory of metals not much attention is usually devoted to orbital phenomena. The atomic orbital 
angular momentum is first of all largely quenched by crystal fields, which split part of the orbital degeneracy. Sub- 
sequently, the electron kinetic energy, usually large, completes the job by eliminating residual degeneracies except at 
special points in the Brillouin zone. 

The scope of this paper is to explore and describe in some detail a possible scenario emerging in the opposite limit. 
This is realized considering "molecular" sites with a well-defined orbital degeneracy, unquenched by crystal fields, and 
with a coherent, but extremely weak, electron hopping energy t between the sites. The specific extra ingredients which 
we consider are: (a) the orbital degeneracy depends on electron occupancy; (b) inter-site electron hopping forces a 
switch of orbital degeneracy for both the initial and final site. 

This kind of model is closely inspired by our previous studies of the physics of Cqq ions, present in molecular 
fuUerides. An isolated Cg^^ ion undergoes a dynamical Jahn- Teller (JT) effect, which, as we found, is affected by an 
electronic Berry phase when n is odd, but not when it is even. As a result, the molecular ground state has L = 1 in 
the former case and L = in the latter as in point (a) above. 

A realistic model for the superconducting compounds AaCgo (A=K,Rb) should take into account the electronic and 
molecular degrees of freedom of Ceo: the three degenerate tig electronic orbitals (LUMO) at each molecule and the 
eight Hg 5-dimensional phonon multiplets which are allowed to couple, by symmetry, with the electronic orbitals. The 
energies of these phonons range from 30 to 200 meV, i.e., comparable to the typical electronic bandwidth W « 0.4 
eV in the solid. In such a situation, a standard Eliashberg approach, heavily based on Migdal's theorem (neglect of 
electron-phonon vertex corrections) and on the smallness of ?iwDcbyc/W^, is hard to justify from first principles, 
even if empirically successful. Q 

We explore here a different route, based on the opposite limit in which the hopping matrix element t is the smallest 
energy scale in the problem. The basic ingredients of the model we will study are the different degeneracies associated 
to molecular states with different electronic occupation, as suggested by a careful treatment of the JT effect of each 
molecule. |H§] While well aware of the limitations of such an approach - the real situations being in the intermediate 
regime t ~ huj - we still believe it illustrates the potentially interesting role played by orbital degeneracies in molecular 
crystals. 

The basic ideas can be illustrated in a model which is much simpler than that of Cgo, i-C, a system of "e (8) E" 
molecules with two electronic orbitals per site c±^a, coupled to a two-dimensional JT active phonon multiplet b± 
of energy Hlo: 

± <r,r'> cr ± r 

+ i^c-ph + (e-e terms) 

r a 

Such a simplified model has been considered in Ref. |^ and studied using a strong coupling {g — *■ oo) approach to the 
electron-phonon molecular Hamiltonian. Consider first the case in which there is no hopping {t = 0). The molecular 
Hamiltonian commutes with an angular momentum-like operator 

Jir) = [blir)b+ir)-bUr)b-{r)] + ^J2^cljr)c+Ar) ^ cljr)c^Ar)] ■ (2) 

a 

One finds that molecules with an odd number of electrons n (n — 3) have an associated electronic Berry phase of 
TT and a four-fold degenerate ground state (j = ±1/2, cr = ±1/2). On the contrary, molecules with an even number 
of electrons (n = 0, 2) have no Berry phase and a non-degenerate singlet ground state {j = 0,S = 0). For g >> 1, 
excited molecular states are separated from the ground state by terms of order TiVl = fiuj / g"^ . ^ 

However, the validity of this analysis is not restricted to a strong coupling regime. For g << \, conventional 
2"''-order perturbation theory in Hc-ph, neglecting retardation effects, gives an effective phonon- mediated electron- 
electron interaction of the type |8| 

F,„, = -S^Y.Y.{cljr)c^ Ar)cl^^,ir)c+^,,{r) + H.c.} . (3) 
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It is simple to verify that the four degenerate n = 1 states Cj_ ^\0) have an energy En^i « whereas the 

lowest n = 2 state is the inter-orbital singlet (c^ ^c[_ — |)|0) with energy i?„=2 ~ —g^huj. ^ Once again, 

(j = ±1/2, a = ±1/2) for the n = 1 states, whereas {j = 0, S = 0) for the lowest n = 2 state. 

Equation ^ shows that the inclusion of the JT effect proves already important in providing a sizable pairing energy 
which acts as a kind of effective negative Hubbard-JJ term, therefore reducing the effect of the local Hubbard repulsion 
between electrons. [Q It is therefore not surprising that, in absence of strong repulsive e-e interactions, this phonon- 
mediated local attraction leads to superconductivity. Manini, Tosatti and Doniach have found that this pairing 
attraction is important at least as long as i < hil. Shelton and Tsvelik have studied, using bosonization in ID, the 
purely electronic model consisting of the hopping term T and the effective non-retarded electron-electron interaction 
He-e in Eq. ^. |^ They too find a strong tendency to superconductivity. However, these negative Hubbard ?7's can 
in practice be easily canceled with Coulomb repulsive Hubbard C/'s. 

We will show here that there are effects, induced by orbital degeneracy, which can lead to superconductivity even 
when repulsive interactions overwin the polaronic attraction. These effects will be found to be more pronounced at 
low carried density. 

As argued above, singly occupied states have an extra degeneracy of orbital origin {j = ±1/2) while the lowest 
doubly occupied state is non-degenerate {j = 0). If the hopping t is much smaller than the lowest local excitation 
energy, we can effectively restrict the Hilbert space at each site by keeping only the ground state for each electronic 
occupation. |^ We are thus led to considering a Hubbard-type model in which each site of the lattice carries an extra 
pseudospin-1 degree of freedom, with the following constraint: a singly-occupied site has, in addition to the usual spin 
degeneracy, a twofold (orbital) degeneracy represented by S^- = ±1, whereas a doubly occupied (or empty) site must 
have = and is therefore non-degenerate, 

Ur^l — > = ±1 , 
nr = 0,2 — > S*; = . (4) 

The electron-pseudospin (EP) working Hamiltonian H is written as follows: 

^ = E E ± H.c.) (S+S;, + H.c.) + Uj^nr^nr^ , (5) 

<rr'> (7 r 

where the S'^'s are spin-1 ladder operators at each site, and the remaining notation is completely standard. It is 
worth stressing that while H, and in particular the hopping term, conserves the constraint in Eq. ^, the model is still 
highly non-trivial even for U — 0. 

Manini, Tosatti and Doniach have introduced and studied this EP model in one-dimension and at half-filling 
By small ring exact diagonalizations and a BCS-like mean field calculation they show that the effect of the correlated 
hopping, even in the absence of any negative-C/ term, is to induce pairing between electrons of opposite spin. 

In this paper we present the results we have obtained for the full phase diagram of the model in ID, which 
is illustrated in Fig. 1. We find that at low densities n << 1 there is a region with dominant superconducting 
correlations surviving repulsive interactions as strong as the bandwidth, W — At. The critical value Uc of U below 
which superconductivity is found to be dominant decreases with increasing density n. For n = 1/2 we find Uc/t ~ 1, 
whereas at half-filling our results (not very accurate, due to marginal umklapp terms) are compatible with Uc/t « 0. 
For U > Uc{n) and out of half-filling (n ^ 1) the system is metallic with dominant 2/cf charge density wave (CDW) 
correlations. At half-filling a charge gap opens up for U > Uc and the system becomes an insulator. A spin-gap 
characterizes the phase-diagram for all densities and for all values of U, even in the metallic regime U > Uc- 



II. SOME DEFINITIONS AND A FEW EXACT RESULTS, MAINLY IN ID 

As we will discuss more in detail in Section 3, only charge degrees of freedom remain gapless in our model. As a 
consequence, the large-distance physics distances is characterized by a quantity Kp - determining the behavior of the 
correlation functions decaying as power-laws - and by the velocity Up of the sound-like collective excitations. ^0 11 



As we will see in Section HI, the relevant correlation functions are density-density {N{x)) and singlet superconductive 
{S{x)), behaving as 

Kp 1 cos {2kpx) cos {Akpx) 



4- 



x'-'^p 
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For general values of U and of the density n, results for Kp can be obtained in a standard way ^ 
of the charge velocity Up, 



from a knowledge 



lim —[EUS 

L— »oo ZTT 



0, A: = 2tt/L) ~ El] 



of the inverse compressibility k 



and of the Drude peak strength Dc [|12|,|13| 



= lim L 



9^ 
dm 



2K„ 



^ lim ^^^^-('^^ 

L — >^oo 2 



-UpKp 



(7) 



(8) 



(9) 



Here El{S = 0, /c = 27r/L) is the energy of the lowest excited singlet state with total momentum k = 27r/L with 
respect to the ground state of energy E^ in a ring of size L, and El{^) is the ground state energy in presence of a 
magnetic flux $. 

In general, we resort to exact diagonalizations of the model on small rings of size L (with L up to 16), and use Eq. |^ 
to estimate Up and Eq. ^ or || to extract Kp. Three regions of the phase diagram are, however, well under control 
using different approaches: the low density limit n —> 0, |15| and the limit of large repulsion or attraction U —> ±oo. 
0,0 We summarize here, for the reader' s convenience, the most important results concerning these three regions. 



A. The low-density limit: solution for two electrons 



Consider the two-electron problem first. A state in the two-particle Hilbert space with total z-component of the 
spin M^Qrp = (for both the electrons and the spin-1 states) can be written as 



= J2 [^P+-{r,r')S+S;, + ^P-+ir,r')S-S+] 



rT 



.t 



|0), 



(10) 



where the vacuum |0) is the state without fermions and with S"^ = at each site. In writing one takes into account 

the two possibilities of associating a S"^ = ±1 spin state to the up and down electrons: Vh is the amplitude for 

having = +1 associated to the t-electron (and S"^ = — 1 to the |-electron), while -0 ^ is the amplitude for other 

possible choice. The Schrodinger equation for ?/'h (r,r') is 



Eijj+^{r,r) = -i ^ [ V-'+-(j' + a, /) + 0+-(r, r' + a) ] + U 5ryilJ+-{r,r') 

a 

-t{J2Sr+a.r') [ V-+ (r, r) (/ , /) ] , 



(11) 



where a denotes a nearest neighbor vector (a = ±1 in ID). A similar equation is obtained for tp ^(r, r') by just 

exchanging -0 (_ and "0-1 everywhere. The last term in Eq. O is the crucial outcome of the extra (orbital) degeneracy 

of singly occupied states, and deserves a few comments. When the two electrons are far enough in the otherwise empty 
lattice, the EP Hamiltonian H simply allows the hopping to a nearest neighbor site of the "composite" object (see 
Fig. 2a) formed by an electron and the associated spin-1 state (first term in Eq. ^l|). In other words, each electron 
retains its orbital index during the hopping process. Things are more subtle when two electrons come to the same site 
r. In such a case (see Fig. 2b), from a doubly occupied site with S*^ = one can reach, upon hopping, two possible 
final states: either each electron keeps its own spin-1 label or the spin-1 labels associated to the two electrons are 
exchanged. It is precisely this second possibility of exchanging spin-1 states that is responsible for the presence of 

tp y in the equation for ip^ and vice- versa (last term in Eq. 11). 

The Schrodinger equation is easily solved (on a square lattice in any D), in momentum space, where it reads: 



[E^2ek]^+-{k) 



(12) 
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Here the case of total momentum P = has been considered, for simpUcity, and is the tight-binding dispersion of the 

free-electron problem (e^ = —2t cos ka)- Introducing the quantity = L^^ (p), and the corresponding 

one for the — -I- amplitude, and imposing self-consistency, we arrive at a simple 2x2 problem 

A; " fc 

(The equation for J y is obtained by exchanging J-| with J y everywhere.) The set of solutions among which the 

ground state is found have = J ^ and their eigenvalues E satisfy the equation 

1 1 _ 2 

k 



It is worth mentioning that in the ordinary Hubbard case, the right hand side of Eq. 14 would simply read 1/U 



A graphical analysis of Eq. 14 readily shows that a bound state solution is present even for U > up to Uc — 'iDt in 
D <2. In D > 3 a finite attractive U is needed to produce a bound state. 

The bound state solution has been worked out analytically in the ID case. The bound state wavefunction 
naturally provides a picture of bound pairs approximately localized on adjacent lattice sites, Remarkably, the 
rather strong attraction responsible for this binding is generated by the kinetic term alone via the presence of the 
additional degrees of freedom. The critical value of the Hubbard repulsion [Uc — 4t in ID] is considerably larger 
than the ground state binding energy at J7 = [Eb/t = (8/-\/3 — 4) ~ 0.618 in ID] showing that this kind of pairing 
mechanism is rather insensitive to the presence of on-site Coulomb repulsion. The same feature is also present in two 
dimensions where Uc = 8t: The enhancement is due to the the larger coordination of the 2D lattice which provides 
an even more efficient delocalization of the electron pair. 

This interpretation of the two-particle ground state, plus the additional finding that there is no phase separation, 
1^ leads to a simple picture, at least for U < Uc, of the low density limit of model (||) both in one and two dimensions: 
the system behaves as a weakly interacting, dilute gas of bosons with an extended core. A superfluid ground state 
must be expected at zero temperature. In one dimension, of course, off-diagonal long range order cannot occur and 
only a long range power-law decay of the density matrix is possible, while in two dimensions a genuine T — Bose 
condensate will form. In terms of the original electrons this implies a standard strong coupling BCS superconducting 
ground state with localized Cooper pairs. A similar scenario has also been proposed in the framework of the one 



dimensional t — J model where bound pairs arc formed at low density |19| and 2 < J/t < 2.95. In that case, however, 
the model is unstable to phase separation, which in fact occurs massively at larger values of J/t, whereas the present 
model shows no such tendency at least in ID. 

From a low-density expansion of the ground-state energy and of the charge velocity an analytic determination of 
the exponent Kp in the limit n —>■ has been given in Ref. [^. It is found that Kp = 2 for U < Uc{n = 0) = At, 
whereas Kp 1/2 for [/ > 4i. |l|] 



B. U ^ +00 limit 



The strong coupling expansion of the EP model has been worked out in Ref. |16|. The resulting Hamiltonian, to 
0{t^ /U), is defined in the subspace of empty or singly occupied sites where electrons are characterized by two internal 
spin- 1/2 degrees of freedom: The usual spin [a) and a pseudospin (r). The latter takes into account the two possible 
states of the original spin-1 rotators allowed by the constraint (^). The resulting Hamiltonian, exactly like in the t — J 
model, contains a kinetic term (with coupling constant —t) and a spin (and pseudospin) dependent contribution which 
scales as t^ /U . In ID, analogously to the Hubbard model case, the total wavefunction factorizes, at [/ ^ oo, for 
all densities: The position of the electrons is determined by a (free) spinless fermion wavefunction while the spin (and 
pseudospin) ordering on the squeezed chain is governed by the effective spin Hamiltonian 



-Jeff 



<ij> 



2a, 



2t, 



(15) 



where the effective coupling constant depends on the electron density n as Jeff = ^n[l — sin(27rn)/27rn]. It turns 
out that, in contrast with the Heisenberg case, the Hamiltonian in Eq. |l^ has a (doubly degenerate) valence bond 
ground state and is characterized by a spin gap of order Jeff. |0 Therefore, we conclude that our model (Eq. |), at 
arbitrary (finite) density has a spin gap for U — » oo. The charge degrees of freedom are, on the contrary, gapless and 
tend to approach the spinless fermion case as U oo. 
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Quantities such as the charge velocity and Drude peak can be obtained analytically for any lattice size L by this 
large-C/ mapping of the charge degrees of freedom to free spinless fermions. The expression for the charge velocity Up 
reads: 

Up{U — +00) = 2t— sin(7r/L) sin(7rn) 2tsin(7m) , (16) 

TT 

whereas the Drude peak strength Dc is given by 

D,(U = +00 = t \ "--^ - sin Trn . 17 
Lsm(7r/L) tt 

As a result, Kp = 1/2 for all densities n 7^ 1 on the boundary line U = +00. Jll| 

C. ?7 — > —00 limit 

The third, rather trivial, region where the physics is under control is that of large on-site attraction. Indeed, in the 
limit U — !■ —00, only doubly occupied and empty sites are allowed. The extra degeneracy of singly occupied sites plays 
therefore no special role, and the model behaves as the corresponding negative U Hubbard model. ||2l| In analyzing 
our data in this limit we will use, for comparison, the numerical results for the Hubbard case obtained from Bcthe 
Ansatz. [|l^ At half-filling, the mapping of the U —00 limit into the antiferromagnetic Heisenberg chain (with 
J = jjjj in the present case |2l]) provides analytic L — 00 results for the Drude peak Dc and the charge velocity Up. 



0,0 We get § 



i?, = 4- + ... = + . (18) 

8 2|C/| ^ ' 

(The dots indicate higher order corrections in t/U .) The corresponding value of Kp at half-filling for large negative 
U is therefore given by Kp ~ 1. ||l^ 

III. THE PHASE DIAGRAM (ID) 

In order to study the model using standard \D Fermi gas techniques, we have to represent the extra spin-1 
degeneracy and the associated constraint in a suitable way. A possible procedure is to first represent the Hamiltonian 
in Eq. ^ in terms of additional hard-core boson operators dra as follows: 

= -i ^ ^ { c\„c,,,d\^,dr',, -f H.c. } + ^ Yy<\<i + "^T"H ] ' (19) 



<rr > cj,o' 

the constraint now being that, at each site, 



< = < . (20) 



Eqs. 19 and ^ give a faithful representation of all the matrix elements of the original model. The idea behind this 
representation is that the extra label 5'^ = -1-1 (or —1) at each singly-occupied site r is represented by a hard-core 
boson d}^^ (or <i|!|). Empty sites have no bosons, whereas the — ^ state of a doubly occupied site is represented by 



two bosons 



In ID, we can transform the hard-core bosons d^a into spinful fermions, by a Wigner- Jordan transformation. |13|| We 
are then led to consider a problem of correlated hopping of two species of fermions in IZ? with a conserved constraint 
(Eq. |2^) . A possible way of treating this problem is to enlarge the Hilbert space by removing the constraint in Eq. 
|20| and to introduce a term T^X^rC^-r ~ '^r)^ with y — > 00 in the Hamiltonian. The correlated hopping term can then 
be regarded as a second-order process obtained from a standard free hopping of the c and d particles as follows: 
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Hv = -i ^ ^{cl^Cr'a + dl^dr'cr + H.C.} + ^^[K^Kl + "rT^rj] 
<rr'> CT r 

+ y^(<-nf)2 + AHv 

r 

P 

^ V ^ Y.^'^r'aCrJra'Cr'a' + Jr^a' dr' } , (21) 
<rr'> a,a' 

where P /V = t, and y is understood to be large. The extra term AHv is included in order to cancel unwanted 
terms generated by second-order perturbation theory in the hopping t. We have thus obtained a reasonably standard 
form of two-chain problem, which we treat in weak-coupling in U and V (assuming continuity with the large V 
regime), by linearizing the band around the two Fermi points at -iikp and introducing right (p = i? = +) and left 
{p = L = —) fermion fields ^ppa<j{x) for every species of fermions (a = c, d, and a =t, [)■ Luckily, the "g-ology" of 
this two-chain problem turns out to be a particular case of the two-chains in presence of transverse hopping treated 
by Fabrizio et al. Q We can therefore directly use the RG equations of Ref. |Q calculated up to third order in the 
couplings (two-loops). Away from half-filling (i.e., when umklapp terms are irrelevant) the result is that the model 
scales to a strong-coupling fixed point and gaps open in all sectors with the exception of the totally symmetric one, 
which remains gapless. Correlation functions can be calculated by bosonizing the strong coupling fixed point. We 
introduce a bosonic representation of the fermion fields 

^paa{x) - ^^e*v/^[Q"-(-)+P*""(-)] , (22) 

where '^aa{x) and Qaa{x) are boson fields such that Iiaa{x) = —'\/Qaa(x) is the momentum conjugate to $acr(a;), and 
a is a short-distance cut-off. In terms of charge/spin and symmetric/antisymmetric combinations of the fields 

Km = ^\^cpi.)±'^M^)]^ (23) 



the bosonized Hamiltonian for the gapless field reads 



1 



n=\j^ dx[upif;(n+)2 + ^(V<i>+)^} + ^[^(^N-N^f+U,KP^], (24) 

whereas the other fields are effectively locked to their strong-coupling values $p = <i>J ~ 0. Here N is the operator 

counting the total number of particles of both species {Nq is its ground state value), and J is the associated current 
operator. ]To[ | We do not calculate Up and K'^ perturbatively, since we are interested in the physical limit V oo. 
We will therefore extract the correct Up and numerically from finite size exact diagonalizations (sec below). 
From Eq. the compressibility relation can be immediately derived. |l0HT^ It reads: 

d^EL _^ Up 

where N = 2Nc is the total number of particles of both species. Similarly, the Drude peak is found to be given by 

Dc = ^^pK ■ (26) 

It is simple to verify that the correlation functions which decay as a power-law at large distances are the 
density-density correlations N{x) = {n^nf)) and the singlet superconductive correlations S{x) = (PJPq) with 
Px = V^cT {x)ipci {x)ipd^ M'ipdi (x). Expressing the fermionic fields in terms of the bosonic ones, and using the Hamil- 
tonian in Eq. e3, we readily find that: 



Kp 1 ^ ^ cos {2kFx) ^ ^ cos {^kpx) 
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We emphasize that all the relations obtained so far, from the compressibilityto the correlation functions, have 
exactly the same form they would take in the negative-?/ Hubbard case ||2^,|2^ if we define Kp = Kfj2. The 
corresponding forms of the correlation functions, compressibility, and Drude peak are then given by Eqs. Ig,^,^- (In 
the compressibility relation, a factor 4 is absorbed by the fact that N — 2Nc, Nc being the number of electrons.) To 
facilitate the comparison with the Hubbard case we will use in the following the notation Kp = K'^ /2, so that the 
latter will be directly comparable to the Hubbard model's Kp. j2^] 

From Eq. we observe that regions with Kp{— K'^ /2) > 1 and Kp < 1 are characterized by dominant correlations 
(i.e., correlations with the slowest decay) of the superconductive and 2kp-CDW type, respectively. We recall that, 
in the low-density limit n — > 0, Kp — i for U > Uc ^ 4:t and Kp = 2 for U < Uc- This defines the left (n = 0) 
border of the phase diagram in Fig. 1. For U < Uc the dominant correlations are superconductive, decaying as l/^/x, 
exactly as in a dilute (n — > 0) hard core boson gas. For U > Uc, the correlation function with the slowest decay is 
the 2kp density-density response, with an inverse square root behavior implying a divergence in the 2kp-CiyW 

susceptibility, except at U = oo where the 2kF response function has vanishing amplitude, like in the Hubbard model, 
and the corresponding singularity is therefore absent. 



A. Numerical results 

We are now ready to discuss the results of the finite-size exact diagonalization data of the EP Hamiltonian (||). 
First we will show that a spin gap - which is expected all over the phase diagram, and not just in the superconducting 
region - is indeed found for all the values of U studied, at n = 1/2 and n = 1. Figs, g and ^ show the finite-size 
gaps between the ground state and the lowest triplet state, AEl{S = 1), as a function of the inverse lattice size 
1/L for both the quarter-filled (Fig. ^ and the half-filled case (Fig. ^) and for several values of U. As expected, 
AEl{S = 1) always extrapolate to a finite non-zero value in the infinite volume limit for all the values of U studied. 
The extrapolated values of the spin gap Ai?oo('S' = 1) (see Fig. ^ - inset for the n = 1 case) decrease with increasing U 
in agreement with the result that, as [/ — > +oo, the spin gap should be proportional to Jeff ~ — sin (27m)/(27r)]. 

Next we turn to the exponent Kp for finite values of the density n. The exact eigenvalues of L-site rings provide 
determinations of the three quantities Up (Eq. |^), Dc (Eq. |9|), and d^EL/dNc (Eq. H). In particular, the magnetic 
flux $ through the ring in the evaluation of the Drude peak (|^) is implemented simply by changing the hopping 
matrix element t — > te'*/^ (a procedure equivalent to twisting the boundary conditions) and calculating the second 
derivative with respect to $ numerically. Moreover, as usual, the second derivative of the energy with respect to 
the number of electrons Nc, appearing in the compressibility (^), is approximated by the finite difference expression 
L[El{Nc + 2) + El{Nc - 2) - 2El{Nc)]/^. 

Figs. ^ and |^ show the finite-size values of the Drude peak Dc and of the charge velocity Up (inset) for the quarter- 
filled [n — 1/2) and half-filled {n — 1) case, for several values of U . The U = +oo values of Dc and Up for n = 1/2 
(see Eqs. ^ and |l^) are shown in Fig. ^ by solid lines. In the quarter-filling case - representative of the < n < 1 
region of the phase diagram - a finite non-zero Drude peak, signal of a metallic ground state, is obtained for any 
U . Correspondingly, the charge velocity Up scales to a finite value for any U . On the contrary, at half filling and for 
positive enough U {n—1, Fig. m the Drude peak vanishes in the thermodynamic limit, indicating insulating behavior. 
This is due to the opening of a gap in the charge sector above some critical value of U close to 0. In Fig. ^ the 
finite-size charge gap at half-filling 

A, = i [E{Nc = L + 2)-E{Nc = L)-U] (28) 

is plotted: Ac scales nicely to as L — > oo for negative U, whereas it extrapolates towards finite nonzero values 
(roughly proportional to U for large U) for positive enough U. The extrapolated values of the charge gap for several 
values of U arc shown in the inset of Fig. ^, together with the positive U Hubbard model values, for comparison. 
Finally, Figs. ^ and ^ show the finite-size estimates of the exponent Kp iQ for n = 1/2 and n = 1, respectively, and 
for several values of U. For the quarter-filled case (Fig. ^ it is instructive to compare Kp with the corresponding 



Hubbard model result, obtained from a numerical solution of the Bethe Ansatz equations. |12| The arrows in panel 
(a) of Fig. P indicate the values of Kp for the Hubbard model at the corresponding value of U. The U/t = —100 data 
confirm that the EP model behaves like the Hubbard model in the U ^ —oo limit. 

The most important information, however, is contained in the [/ = data of Fig. |[ the exponent Kp atU — O scales 
to a value somewhat larger than the critical value 1, which is also the value of Kp for the non- interacting Hubbard 
model. This indicates (see Eq. with Kp — K^ /2) that our model has dominant superconducting correlations even 
in absence of attractive interactions. Panel (b) of Fig. shows that the superconducting state survives, for n = 1/2, 
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up to a positive U /t of order 1. Fig. ^ shows the values of Kp at half-filhng, together with the Hubbard model results 
(solid squares). At half-filling, umklapp processes are marginal even in the gapless regime, introducing logarithmic 
corrections to the size scaling. |2^] The data for the Hubbard model show clearly the consequence of such a non-trivial 
size scaling: an extrapolation based on the data at small sizes gives sizable errors (of the order of 10%). In view of 
this non-trivial size scaling, it is hard to make definite statements on the critical value of U above which a charge 
gap opens and the system starts to be an insulator. Nevertheless, the U — data for Kp (panel (b) of Fig. ^ seem 
to be definitely below the value of 1, where umklapp processes should become relevant. We are therefore inclined to 
believe that a charge gap should be present even at J7 = (see also the comparison with the Hubbard results of Fig. 
1^, although the U = Q Drude peak results of Fig. || might suggest a conducting behavior. 

The phase diagram in Fig. |l| summarizes our results on the model in Eq. (^) in Z? = 1. The most remarkable feature 
is in the low density region n ^ 1, where superconducting correlations are found to be dominant up to positive values 
of the on-site coulomb repulsion U comparable to the electronic bandwidth. We believe that the physics at small 
densities is robust. The bound states discussed in Section [I A should lead to Bose condensation of strongly localizes 
Cooper pairs in 2D. In ZD no bound state is found but a superconducting instability should still occur in the presence 
of a Fermi surface. On the other hand, other features of the model, like the presence of a spin gap over all the phase 
diagram - including the positive U metallic regime - are probably confined to the low dimensionality ID case. 



IV. CONCLUSIONS 



We have studied the \D phase diagram of a model of electrons which hop weakly between molecular sites. Due to 
the Berry phase mechanism, each molecular site possesses an unquenched orbital degeneracy, only when its electron 
occupancy is odd. Upon electron hopping, the orbital degeneracies must correspondingly switch, as the occupancies 
of the initial and final states change. By adding adjustable on-site Hubbard U term, we can effectively remove from 
the problem all polaronic effects, and we can thus study the effect of orbital degeneracy, and of its hopping-induced 
switch, in isolation. The working model (|^) is therefore Hubbard-like, but embodies the orbital degeneracy part in 
the form of a pseudospin-1 on-site variable. We have studied its zero-temperature phase diagram as a function of 
electron filling and of Hubbard U . 

Through a series of mappings - all of which are exact for what concerns the asymptotic behavior of correlations 
in ID - we find that the problem is equivalent to a particular variant of two coupled 11? chains, whose behavior is 
generally well studied . This mapping predicts that there will generally be no charge gap (except at half filling for 
U > Uc) while there will always be a spin gap, for all n > and all U. The zero charge gap is related to either sliding 
CDW conductivity or superconductivity, and we have distinguished these two different regions in the phase diagram. 
A quantitative study of the phase diagram, conducted by exact diagonalization and finite size scaling, has led to the 
overall picture of Fig. [|. 

The most interesting fact in the phase diagram is the dominance of superconductive behavior extending even into 
the (moderately) repulsive-?/ region. Pairing is associated with new orbital correlations - as foreshadowed by an exact 
2-electron solution - which permit a lowering of kinetic energy through an enhanced nearest neighbor pair scattering. 
This feature represents an important novelty brought about by orbital degeneracy, and should reasonably carry over 
to more complex and realistic situations in higher dimensions. This pairing mechanism is not easily destroyed by a 
repulsive U, it is more effective at low carrier density, and is apparently immune from the polaron self-trapping, which 
depresses in strongly coupled electron-phonon systems. 

Another new feature found in this model in ID is the presence of an overall spin gap, including the nonsupercon- 
ductive region of the phase diagram, where sliding CDW 2kp correlations dominate. It is at this stage not clear in 
what way this behavior (which appears rather specific of ID) will change in higher dimensions. 
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FIG. 1. The ID phase diagram for the EP model (Eq. ^ discussed in the present paper. For details see Sect. |l|. 
FIG. 2. The hopping mechanism of an up-electron starting from a singly (a) or doubly (b) occupied site. 

FIG. 3. The spin-gap, defined as the difference in energy between the lowest triplet state and the ground state (singlet), for 
n = 1/2 and U = 0,2t and 4t, as a function of the inverse lattice size 1/L 

FIG. 4. The spin-gap for n — 1 and U/t = 0,1,4,10 and 100. Inset: the extrapolated values of the spin-gap in the 
thermodynamic limit. The extrapolations are obtained through weighted quadratic polynomial fit to finite size data. The 
asymptotic 1/U behavior is apparent at the right side of the plot. 

FIG. 5. The Drude peak Dc and the charge velocity Up (inset) for n — 1/2 (qua rter- filling). The lines for the Drude peak 



are fit to the finite-size data. The large-17 behavior (solid line) is discussed in Sect. II B 



FIG. 6. The Drude peak and the charge velocity Up (inset) at half-filling (n = 1). The data for U = 2t and lOt show 
clear indication of Dc vanishing in the thermodynamic limit, while for negative U , Dc remains finite. The data for U = Q seem 
also to converge to a nonzero value, but there are evidences (see Sect. |tl|) in favour of insulating behavior {Dc ~ 0). 
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FIG. 7. The charge gap at half filhng, obtained as the difFerence in energy between the lowest Nc = L -\- 2 state and the 
Nc = L ground state (Eq. |2^). The inset reports the extrapolations obtained through weighted quadratic polynomial fit to 
finite-size data, together with the Hubbard model values, for comparison. The charge gap seems to vanish at some point close 
to ?7 = 0. 

FIG. 8. The exponent Kp for the EP model at quarter filling obtained (Eq. ||) as Kp = irDc/up. Panel (a): U/t — —100, —4 
and 0; arrows mark the asymptotic values of Kp for the Hubbard model obtained by Bethe Ansatz. Panel (b): U = t and 4t; 
squares mark the scaling of Kp for the Hubbard model at the corresponding values of U, for comparison. 



FIG. 9. The exponent Kp for the EP model at half filling obtained (Eqs. | and ^ as Kp ^ n (^Dc/2L^^^^j . Square 

dots represent the scaling of the same quantity for the Hubbard model at the corresponding values of U. Notice the non-trivial 
size scaling, due to marginal umklapp terms. Panel (a): U — —20t and ~2t. Panel (b): U — and t. 
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